arXiv:1507.03094v2 [hep-lat] 28 Jan 2016 


Prospects for a lattice computation of rare kaon decay amplitudes 

K —>■ decays 

N.H. Christ/ X. Feng/ A. Portelli/ and C.T. Sachrajda^ 

(RBC and UKQCD collaborations) 

^Physics Department, Columbia University, New York, NY 10027, USA 
‘^School of Physics and Astronomy, University of Southampton, Southampton S017 IBJ, UK 

(Dated: January 29, 2016) 


Abstract 

The rare kaon decays K —)• and K —)■ ttuV are flavor changing neutral current (FCNC) 

processes and hence promising channels with which to probe the limits of the standard model and 
to look for signs of new physics. In this paper we demonstrate the feasibility of lattice calculations 
of A —7- decay amplitudes for which long-distance contributions are very significant. We 

show that the dominant finite-volume corrections (those decreasing as powers of the volume) are 
negligibly small and that, in the four-flavor theory, no new ultraviolet divergences appear as the 
electromagnetic current J and the effective weak Hamiltonian Hy/ approach each other. In addition, 
we demonstrate that one can remove the unphysical terms which grow exponentially with the range 
of the integration over the time separation between J and Hw- We will now proceed to exploratory 
numerical studies with the aim of motivating further experimental measurements of these decays. 
Our work extends the earlier study by Isidori, Turchetti and Martinelli [1] which focussed largely on 
the renormalization of ultraviolet divergences. In a companion paper [2] we discuss the evaluation 
of the long-distance contributions to A —?■ ttuV decays; these contributions are expected to be at 
the level of a few percent for A"*" decays. 
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I. INTRODUCTION 


The rare kaon decays K —)■ and K —)■ ttvV are flavor changing nentral cnrrent 

(FCNC) processes which arise in the Standard Model throngh W-W and 7 /Z-exchange 
diagrams, containing up, charm and top quarks in the loop. As a second-order electroweak 
interaction, the SM contributions are highly suppressed in FCNC processes, leaving the 
rare kaon decays as ideal probes for the observation of New Physics effects. Additionally, 
these decays can be used to determine SM parameters such as Vu and Vts, to search for 
CP violating effects and to test the low-energy structure of QCD as described within the 
framework of chiral perturbation theory (ChPT). In this paper we discuss the feasibility of 
computing K —)■ decay amplitudes in lattice simulations; the corresponding study for 

K —)■ TTuV decays will be presented in a forthcoming companion paper [2]. 

The first observation of 41 —)■ decays was made at the CERN PS accelerator 

in 1975 [3]. After a long series of experiments spanning 40 years, NA48/2 at the CERN SPS 
accelerator has observed a sample of 7253 —)■ vr^e^e” decays [4] and a sample of 3120 

—)■ decays [5]. These precision measurements give important information on the 

low-energy structure of the weak interaction and provide important tests of ChPT. Though 
expected to be difficult, the first observations of the decays Ks —t 7 r°e+e“ (7 events) [ 6 ] 
and Ks —t (6 events) [7] are reported by the NA48/1 experiment at the CERN 

SPS accelerator. These Ks decays are important in isolating the contribution of direct 
CP violation in the decay Kl ^ For these interesting CP-violating Kl decays, 

upper bounds are set for Br(iCi —)■ 7 r°e+e“) < 2.8 x 10 “^° [ 8 ] and Br(iFi —)■ 7 r°p+/i“) < 
3.8 X 10-1° [9]. 

On the theoretical side, much work has been done to understand and evaluate both the 
short- and long-distance contributions to rare kaon decays. Some useful reviews can be found 
in [10-15]. The CP-conserving decays K^ —)■ and Ks —t are dominated 

by long-distance hadronic effects induced through the one-photon exchange amplitude. So 
far the relevant decay amplitudes are studied and parameterized within the framework of 
ChPT [16, 17]. A challenge, but also an opportunity, for the lattice QCD community 
is to compute the decay amplitudes reliably, as well as determining the necessary low- 
energy constants used in ChPT. We will explain in the following section that the significant 
interference between the direct and indirect CP violating components of the decay K^ —)■ 
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(see Refs. [18, 19]) implies that lattice QCD results for Ks decays can be used 
to evaluate the CP-violating contributions to Kl decays. In Ref. [1] it had already been 
proposed to calculate the long-distance contributions to rare kaon decays using lattice QCD. 
Our work builds on Ref. [1] and leads us to conclude that such computations are feasible 
with present understanding and recent theoretical and technical advances. 

The remainder of this paper is organized as follows. In the next section we review the 
phenomenological background for K —)■ decays with an emphasis on the long-distance 

contributions. The procedure to obtain the amplitudes from lattice simulations is presented 
in Sec. III. Finite-volume effects and renormalization are briefly discussed in Sec. IV and V 
respectively. Finally we present our conclusions in Sec. VI. 


II. PHENOMENOLOGICAL BACKGROUND 

In A —>■ decays, the loop function associated with the y-exchange diagrams has 

a logarithmic dependence on the masses of the quarks entering in the FCNC process (this 
is known as logarithmic or soft Glashow, Iliopoulos, Maiani (GIM) mechanism). The un¬ 
suppressed sensitivity to the light-quark mass is a signal of long-distance dominance in the 
GP-conserving —>■ and Ks —)■ 'n'^K£~ decays. The short-distance contribution 

to the amplitude from Z-exchange and W-W diagrams also exists, but is much smaller than 
the long-distance part induced by the y-exchange diagrams and does not play an important 
role in the total branching ratio. This logarithmic GIM mechanism does not apply to direct 
GP violation in Kl —)■ decays since ImA„ = 0 where Xg = V*gVqd. As a result, the 

direct GP-violating contribution is short-distance dominated. 

Gonsidering only the dominant one-photon exchange contribution, the amplitude A, {i = 
•F, S) for K^ and Ks decays can be written in terms of an electromagnetic transition form 
factor V(z) [17, 20] 

A = - Vi(z)(k + p)>^ ue(p_)7f^ye(p+), ( 1 ) 

where and q = k — p. Here we follow the notation used in Ref. [15]; k, p and p± 

indicate the momenta of the K, tt and respectively. For both K^ and Ks decays, the 
form factor Vi{z) has been analyzed in GhPT [17] and has been parameterized in the form 

V,{z) = a, + hz + Vr{z). ( 2 ) 
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Figure 1. One-loop contribution in CliPT to the decays Ks —t and K~^ —)• . The 

dashed, circled and wiggly lines represent the pions, kaon and photon respectively. There is a 
branch cut when > 4:M^. 

Here Vi(z) is analytic in the complex 2 ;-plane, with a branch cnt starting from 4r^, where 

= M-^/Mk- As shown in Fig. 1, at low energies the tttt intermediate state is expected to 
play the dominate role. Thus is introduced to take account of the 7 * —)■ vrvr effects. 

The contribution of excited intermediate states is not given explicitly and may be accounted 
for by the assumed polynomial correction Oj -|- biZ. A detailed expression for V^'^^(z) is given 
in [17]. As a standard twice-subtracted dispersion relation, V^'^'^(z) satishes = 0 . 

Therefore, the inclusion of V^^^(z) does not affect a*, which is the form factor Vi(z) at zero 
momentum transfer. The a* and bi, can be determined using the experimental data using the 
dilepton invariant mass spectra as inputs. The parameterization (2) provides a successful 
description of the —)■ data but shows large corrections beyond leading order 

in ChPT [17]. A lattice QCD calculation can help to understand the origin of these large 
corrections. 

For Ks decays, only a few events have been observed in experiments. The dilepton 
invariant-mass spectra are therefore unavailable. Assuming vector meson dominance, the 
authors of [17] used the branching ratios to determine the parameter la^l and obtain |a 5 | = 
1.06lg2i the electron and |a 5 | = l.Sd^gg® for the muon. There are two drawbacks of 
using the branching ratios: they do not give information about the explicit z-dependence 
and only the modulus of 05 can be determined. It would be very useful if lattice QCD 
calculations could determine the sign of as and also provide a test for the vector meson 
dominance assumption. 

Ks —)■ decays are interesting for precision studies of CP violation. The relevant 

decay amplitudes receive three major contributions [18, 19]: 

1. a short-distance dominated direct CP-violating term, see Fig. 2(al,a2), 
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Figure 2. The three major contributions to Kl —>■ decays: (al,a2) short-distance dominated 

penguin and box diagrams, (b) long-distance dominated indirect CP violating contribution occur¬ 
ring through mixing, (c) long-distance dominated CP conserving component involving 

two-photon exchange. 


2. a long-distance dominated indirect CP-violating contribution from the decay of the 
CP-even component of Kl (iCi ~ Ks), see Fig. 2(b), 


3. a CP-conserving component which proceeds through two-photon exchanges, see 
Fig. 2(c). 


These three contributions are of comparable size [18]. Here we mainly focus on the CP violat¬ 
ing effects. The total CP-violating contributions to the Kl branching ratios are summarized 
in [15] as 


Br(iC£, —)■ )cpv = 10 x 

Br(iCj^ —)■ )cpv = 10 x 


15.7|as|^ ±6.2|as 


/Im At 
V 10-4 


+ 2.4 


/Im At 
V 10-4 


2 


3.7|as|2±1.6|a5 


/ Im At 

Vio-4 


+ 1.0 



( 3 ) 

( 4 ) 


where the lagp-term comes from the indirect CP-violating component of the amplitude, the 
(ImAt)^-term comes from the direct CP-violating component and the |a 5 |(Im At)-terms are 
the interference between the indirect and direct components. The ± symbol indicates that 
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the sign of as is unknown. |a 5 | are quantities which are expected to be of 0(1) and the CKM 
matrix element takes the value ImAt/10“'^ ~ 1.34. A change in the sign of as can cause a 
large difference in the predicted branching ratios for both —)■ 'K^e^e~ and Ki ^ vr^/i■*■/!“ 
decays. Once lattice QCD has determined the sign of as, this large uncertainty will be 
clarihed. 

In this paper we do not address the evaluation of the CP-conserving (CPC) component of 
the iCn —amplitude given by the two-photon exchange diagram in Fig. 2(c). The he- 
licity suppression of the angular momentum J = 0 channel leads to a negligible contribution 
(of 0(10“^^) to the branching ratio) for the electron mode, but a comparable contribution 
to the CP-violating one for the muon mode. For example, using a phenomenological study 
based on ChPT, the authors of Ref. [18] estimate 

Bi{Kl ^ 7rV+/U")cPC = (5.2 ± 1.6) lO"^' . (5) 

The J = 2 contribution on the other hand is expected to be negligible for the muon chan¬ 
nel, but may be of 0(10“^^) for the electron channel so that the decay Kl —>■ T:^e^e~ is 
predominantly CP-violating [15]. Nevertheless, in due course after we manage to compute 
the Kl ^ —)■ contribution to the decay amplitude, it will be an interesting 

challenge to compute the two-photon exchange contribution. 

Since the phenomenology of rare kaon decays has, up to now, been based on ChPT much 
of the discussion of this section has been in this context. We stress however, that the goal 
of lattice computations reaches beyond the evaluation of the low energy constants. The 
amplitudes will be calculated from hrst principles at several values of and the results 
themselves can then be used in future phenomenological studies. 


III. EVALUATION OF THE AMPLITUDES IN EUCLIDEAN SPACE-TIME 

In this section we discuss the evaluation of rare kaon decay amplitudes using lattice 
computations of Euclidean correlation functions. We start however, with the dehnition and 
a discussion of the amplitudes in Minkowski space. 
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A. Definition of the amplitude 


The long distance part of the K'^ —)• 7r~^£~^£ and —)■ n^£~^£ decay amplitudes is given 

by |17, 20|: 

= /d'l' (ir^(p)l T14(0)i/„,(x)| |/e(k)) . (6) 

The external states are on their mass shells and we define q = k — p. The index j = +, 0 
labels the charge of the mesons, and Hw is the effective weak Hamiltonian density defined 

by |1|: 

HM:^) = ^ v;,v^ |C,(Or - «;) + - Q ‘,)], (7) 

where the Ci are the Wilson coefficients in a chosen renormalization scheme. Ql 2 are the 
following current-current local operators: 

<5? = {Sa7ida){qb7^'"qb) and Ql = {sa^iliqaWbl^^db) , (8) 

where a and b are summed color indices and 7^ = 7/^(1 — 75 )- We envisage working in the 
four-flavor theory and exploiting the GIM cancellation of ultraviolet divergences as explained 
in Sec. V. The electromagnetic current is the standard flavor-diagonal vector current: 

4 = ^( 2 r;-y' + 2 V'-;-\/‘). ( 9 ) 

Electromagnetic gauge invariance implies that each non-local matrix element can be writ¬ 
ten in terms of a single invariant form factor: 

f d-i (7r'(p)| T|J„(0)(Q“(i) - Qt(x))] |Ab(k)) = [g^k + p), - (Ml - M^q,] 

(10) 

and the non-perturbative QCD effects are contained in the form factors wl (g^). In the phe¬ 
nomenological studies described in Sec. II, the wl (g^) are written in terms of a parametriza- 
tion influenced by ChPT. Note that a consequence of (10) is that the matrix elements vanish 
when p = k = 0, i.e. when the invariant mass of the lepton pair is the largest which is 
kinematically allowed, g^ = g^ax- Thus to obtain non-zero matrix elements at least one of 
the mesons must have a non-zero three-momentum. 

Inserting a complete set of states in each of the two possible time-orderings in (6), one 


7 



obtains the following spectral representation for the amplitude: 



• .p P(^) (^^'(p)l MO) \E,k) {E,k\ Hw{0) \K^{k)) 

Jo 2E EK{k)-E + te 

.p (7r^'(p)l HwjO) \E,p) {E,p\ J,{0) \K^{k)) 
Jo 2E E-E^{p)+ie 


( 11 ) 


where p and ps are the associated spectral densities. For what follows it is important to 
notice that p {ps) selects only states with strangeness 5 = 0 (F* = 1). 


B. Euclidean correlators 

For the remainder of the paper we assume that the vector current is the conserved one as 
given by Noether’s theorem and which depends on the chosen lattice discretization of QCD. 
The rare kaon decay amplitudes will be determined by computing Euclidean correlation 
functions and we now turn to a discussion of these. Throughout this section we consider the 
time dimension to be inhnite or very large compared to the time separations of the inserted 
operators. 


1. 2-point functions 


To obtain the energy of single-meson states we consider the following 2-point function; 


r®(«,p) = (^p(«,p)^p(o,p)>, 


( 12 ) 


where (j)p{t,p) is an annihilation operator for a pseudoscalar meson P with spatial momen¬ 
tum p at time t. For t 3> 0 (and t <^T/2, where T is the temporal extent of the lattice), 
Fp^(t,p) has the following behavior; 



-Ep{p)t 

2Epip) 


(13) 


with Zp = (0| 0p(O,O) |P(Ep(p),p)) and Pp(p) = \/Mp -|- p^. Depending on the choice of 
interpolating operator <pp, Zp may (and in general will) depend on the momentum p but 
we do not exhibit this dependence explicitly here and in the following. 
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Figure 3. Diagrams contributing to the 3-point function {tH,p) defined in Eq. (14). The two 
black circles represent the currents in the four-quark operators Qf 2 defined in (8). i denotes a light 
(u or d) quark propagator. The different topologies contain the operators Q\ or depending on 
whether the initial state is a charged or neutral kaon. For example, when the initial state is , the 
W and S topologies contain an insertion of while the C and E topologies contain an insertion 
of Ql 

2. 3-point functions 

To obtain matrix elements of the effective weak Hamiltonian, we define the following 
3-point function: 

(%,p) = j d^x (14) 

with 0 < tfj < and on a discrete lattice the integral over x is replaced by the corresponding 

(3) 

sum. The 4 Wick contractions necessary to compute are illustrated in Fig. 3. For 
(3) 

0 ^tn ^ t-K, Th has the following behavior: 

4-C/7r(p)F/A'(p) 

with ^h{p) = {7r(p)| Hw{0) |-fi"(p))- We also dehne the 3-point function of the electromag¬ 
netic current: 

= y d^xe“*'i'^ ((/)p,(f,p)J/,(fj,x)0^(O,k)) , (16) 
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where P denotes the pseudoscalar meson (P = vr or K) and j its charge. This correlation 
function has the following asymptotic behavior for f S> tj 3> 0: 

(P,k) 


.(3)PJ 








4Ppj(p)Ppi(k) 

where .4'_^’(p,k) = (P’(Epj(p), p)l Jp(0) |P.'(Bpj(k), k)) (note that .4'^°(p,p) = 0), 


(17) 


3. A-point functions 


In order to compute the amplitude (6), we define the following unintegrated 4-point 
correlation function; 

rif^^(tHpj,k,p) = y d^x J dVe“*'^'''(0^i(^^,p)T[J^(tj,x)PH/(trr,y)](/>x. (O,k)) , (18) 

where 0 < tj, tn < tn- As explained in the next section, the rare kaon decay amplitudes are 
obtained by integrating tj, k, p) over tn and tj (or by exploiting time translation 

symmetry and integrating over their difference). 

We now perform the quark Wick contractions in (18) to generate the diagrams which 
need to be evaluated. Assuming isospin symmetry in the quark masses, = m^, 20 types 
of diagrams have to be computed for the charged correlator and 2 additional ones are needed 
for the neutral correlator. We organize these diagrams in 5 classes, which are presented in 
Figs. 7-11. It is convenient to define the factor 


Z'K.Txit-jp ^ k, p) 




3 E-r: (p)^?! 


(19) 


4£,(p)Bp-(k) 

which represents the propagation of the external pseudoscalar mesons in tj, k, p). 

This factor does not contribute to the rare kaon decay amplitude and we choose to define the 
normalized unintegrated correlator jZKp. The decay amplitudes are obtained by 

integrating over tn and tj as explained in the following subsection. We note however, 
that if the times are sufficiently separated for to be dominated by single particle 

intermediate states, then one has: 




rlf^^(^H,o,k,p) = < 


2E^{k) 


if 0 tn to 


^h{p)^j (P, k) ^-EKik)tj^-EK{p)(tH-tj)^EAp)tH if ^ ^ _ 

2Ek{p) 

( 20 ) 
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C. Extracting the rare kaon decay amplitude 


In order to obtain the amplitude (6), we need to integrate the 4-point correlator 


44) i 


{tH,tj,k,p), defining the integrated correlator by: 


P(T„T,,k,p) = e 




rtj+Tb 




P) 


( 21 ) 


^J~Ta 


where Ta,Tb > 0. For T^, Tb such that 0 t j — < t j -|- T;, ^ this integrated correlator 

has the following spectral representation: 


P(T„,T,,k,p) = 


p{E) (7r^(p)| J,(0) \E,k) {E,k\ Hw{0) \K^{k)) 


'0 


2E 


EKik) - E 


(1 _ g[i?x(k)-E]r,^ 


r'+oo 


+ 


dE 


ps{E) {Ti^{p)\Hw{d) \E,p) {E,p\ J,{0) \K^{k)) 


(1 


,-[E-E, 


(P)m)^ (22) 


2E E - E^{p) 

where, as in Minkowski space in Eq. (11), p and ps are the spectral densities of non-strange 
and strangeness S = 1 states respectively. In finite-volume, we write the spectral densities 
as p{E) = ^^(2E„)5(E — En) (and similarly for ps) so that the integrals reduce to sums 
over the finite-volume states n. The exponential factor in Eq. (21) is introduced to cancel 
the tj dependence in (22). To recover the Minkowski amplitude (6), one needs to consider 
the Ta,Tb —)■ -|-oo limit of P^(Ta, Tb, k,p). Since p selects states with strangeness S' = 0, the 
contribution of the states with E < Ex(k) diverge exponentially as —)■ -|-cx) in the first 

integral of (22). This is a standard feature in the evaluation of long-distance contributions 
(see e.g. [21, 22] for a detailed discussion in the context of the K^-Ks mass difference). 
These contributions must be subtracted in order to extract the rare kaon matrix element 
from Euclidean correlation functions. Note that there are no exponentially growing terms 
with Tb, since all the strange states with momentum p have energies larger than E.„-{p). 
We define the subtracted, integrated correlator /^(T^, Tf,, k, p) by subtracting all the terms 
which grow exponentially with from the right-hand side of (22). With this definition we 
can write the rare kaon decay amplitude as follows: 


■<^(9^) = fjTa,Tb,k,p). 




(23) 


How does one compute the subtracted quantity /^(Tq, T;,, k, p) in practice? For physical 
values of the quark masses, the only intermediate states that can generate an exponentially 
growing term in (22) are ones consisting of 1, 2 or 3 pions, the vacuum intermediate state 
being forbidden by parity conservation. We will now discuss each of these cases, providing 
two different approaches for the treatment of the dominant single-pion contribution. 
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1. Removal of the single-pion divergence, first method 


The unphysical divergent term in (22) coming from the single pion intermediate state is 
given by: 

1 


L>^(r„,k,p) = 


,[i?x{k)-E^(k)]r, 


(24) 


2EM EM - EM 

where the notation of Sec. IIIB has been used. Because the pion is a stable state in QCD, 
T)^(Ta,k, p) can be entirely determined by fitting the asymptotic time behavior of the 2 and 
3-point functions as described in Sec. Ill B. 


2. Removal of the single-pion divergence, second method 

We propose a second method to remove the single-pion divergence which is more “auto¬ 
matic” than the previous approach. This method is analogous to the procedure used in the 
calculation of the K^-Ks mass difference in [22]. It is based on an additive shift of the 
weak Hamiltonian: 

k) = Hw{x) -h C5(k)s(x)d(a;), (25) 

where Cs(k) is chosen such that: 

(vr^(k)|W(,(0 ,k)|W^(k))=0 . (26) 

By replacing Hw with H[y the divergent contribution from the single-pion state (24) is 
cancelled. We now show that the transformation (25) does not affect the rare kaon decay 
amplitude. The scalar density appearing in (25) can be written as a total divergence using 
the following vector Ward identity (which is satisfied exactly on the lattice): 

i{ms - md)sd = ) (27) 

where is the sd flavor non-diagonal vector Noether (conserved) current. The relevant 
matrix elements of the scalar density can now be written as : 

(7r^(p)| s(x)d(a;) |E,p) = i- — {p)\Vq^{x) \E,p) (28) 

rris - rrid 

{E, k| Md{x) \K^{k)) = i ^K(k)-E \K^{k)) . (29) 

rris - rud 
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Figure 4. Two-pion intermediate state contribution to the rare kaon decay amplitude. The dotted 
and dashed lines represent respectively the kaon and pion propagators. 

Using (22) and (23) we find that the total contribution of Cs(k)s(i to the amplitude 
is proportional to: 

j d^xe"*‘i'^(7r^(p)| [J/,(tj,x),Q"‘^] |ii:^(k)) = 0 (30) 

because of the vanishing commutator between the flavor-diagonal current and the flavor 
non-diagonal vector charge = J d^y VQ‘^{tH, y)- Thus the physical amplitude is invariant 
under the transformation in Eq. (25). This property is independent of the value of C 5 (k) (and 
thus from the tuning condition (26)). 



3. Removal of the two-pion divergence 


In principle, a two pion intermediate state can contribute to a rare kaon decay through 
the process illustrated in Fig. 4. The matrix elements of vector and axial currents between 
a single-pion and a two-pion state have the following form factor decomposition: 


(vr(pi)l 14 | 7 r(p 2 ) 7 r(p 3 )) = e^upaPlP 2 P 3 F{s,t,u) (31) 


where s = (pi -k ^ 2 )^, t = {pi - p^Y and u = {p 2 - PsY- 
We now show that the vector current does not contribute. Indeed, in Fig. 4 the tttt —>■ vrq* 
vertex gives the following factor: 


^fiupaP k 


dH 




(32) 


(2ir)‘ {P + Ag)l{k - ep + M2] ■ 

Because of 0(4) invariance the integral in (32) can only be a linear combination of p^ and 
which gives a vanishing contribution once contracted with the Levi-Civita symbol. 

In the lattice theory, the cubic symmetry is sufficient for the integral (or the corresponding 
sum in a finite volume) to be a vector, but with corrections which vanish as the lattice spacing 
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a —)■ 0. At finite lattice spacing however, there is a non-zero two-pion contribution from 
lattice artifacts. For example, since the four-component quantity (/c 2 )^, (^ 4 )^} 

transforms as the same four-dimensional irreducible representation of the cubic group as 
A;, one can imagine terms of the form be present. These terms will be 

amplified by the growing exponential factor in (22) and will need to be considered in the 
analysis. By studying the behavior with a? and Ta we anticipate being able to confirm 
our expectation that these effects are very small. For example, in our study of AM^, the 
Ki-Ks mass difference [21, 22], with an inverse lattice spacing of 1.73 GeV and a pion mass 
of 330 MeV, we find that the on-shell two-pion contributions are just a few per-cent and 
the artifacts are of 0(3%) of these. Assuming similar factors here, the exponential factor 
g[Ex(k)-£;]Ta (^ 22 ) would be insufficient for practical values of Ta to make the two-pion 

contribution significant until the calculations reach sub-percent precision. 


4- Removal of the three-pion divergence 

Contributions containing three-pion intermediate states are generated by diagrams such 
as those in Fig. 5. By comparing the measured widths of Ks —)■ tttt decays to those of 
Ks,+ —t vrTTTr decays we estimate the relative phase-space suppression to be a factor of 
0(1/500) or smaller. Moreover, as explained above, we already expect the on-shell two-pion 
contribution to be very small (of order a few percent) and so we anticipate that the on-shell 
three-pion contribution is negligibly small. 

For the the diagram in Fig. 5(a) the contribution to the growing exponential in (22) 
can be avoided completely by restricting the calculations to < 4M^, thus cutting out 
a small region of phase-space. This still allows us to determine the amplitudes in most of 
the range and to compare lattice results with ChPT-based phenomenological models and 
data where this is available. Although it is the diagram in Fig. 5(a) which is dominant in 
phenomenological analyses based on ChPT [17], the imaginary part, corresponding to the 
three-pion intermediate state, is neglected.The exponentially growing terms from diagrams 
such as Fig. 5(b) cannot be avoided in this way and we rely on the phase-space suppression 
described above. (Moreover, much of the exploratory work necessary to develop control of 
the different aspects of the procedure described in this paper, will be performed at heavier 
u and d quark masses, with the mass of the kaon below the three-pion threshold.) 
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Figure 5. Examples of contributions from a three-pion intermediate state to rare kaon decays. 

In the relatively distant future, if the precision required by experimental measurements 
and achievable in lattice computations is sufficiently high one can imagine, in principle 
at least, reconstructing the contributions from three-pion intermediate states explicitly as 
proposed in Sec. Ill Cl and removing the associated exponential divergences in (22). Com¬ 
puting the relevant matrix elements is very challenging however, and the recently developed 
theory of trihadron states on a torus [23] is, so far at least, purely theoretical and has not 
been used in a practical lattice calculation. K —)■ tttttt matrix elements have also not been 
computed to date. 


IV. FINITE VOLUME EFFECTS 

In this section we briefly discuss issues concerning the finite-volume corrections arising 
from the evaluation of rare kaon decay amplitudes computed in a finite, periodic hypercubic 
lattice. The creation of on-shell intermediate multi-particle states will generate finite-volume 
corrections which decrease only as powers of the volume and not only exponentially. Thus in 
addition to generating exponentially growing terms in (22) (see the discussion in Sec. IIIC), 
power corrections are present if there are multi-hadron intermediate states with a smaller 
energy than Mk- We therefore have to identify the potential on-shell intermediate multi¬ 
particle states which can be created in rare kaons decays. This is similar to what was done in 
Sec. Ill C for the subtraction of unphysical divergences in the Euclidean amplitude. We have 
already shown in Sec. Ill C 3 that with a vector current there are no two-pion intermediate 
states and that for < 4M^ there is no power correction from diagrams such Fig. 5(a). 
The arguments given in the preceding section that the remaining contributions from three- 
pion intermediate states are negligibly small applies here as well. Thus, at the levels of 
precision likely to be achievable in the near future, we do not have to correct for power-like 
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finite-volume effects. 

There has been considerable work recently devoted to understanding the hnite-volume 
corrections in three-hadron intermediate states [23]. If and when the precision of lattice 
computations of rare kaon decays amplitudes reaches the precision requiring the control of 
these effects, then it is to be hoped that the theoretical understanding provided in [23] can 
be developed into a practical technique generalizing the use of the Lellouch-Liischer factor 
in it' —)■ TTTT decay amplitudes [24]. 

V. RENORMALIZATION 

The ultraviolet divergences which appear in the evaluation of the matrix elements of 
the form f d^x (/| T[Oi( 0 )O 2 (a:)] |i), where O 12 are local composite operators and |i), |/) 
represent the initial and hnal states, may come from two sources. Firstly O 12 themselves 
generally require renormalization and secondly additional divergences may appear as the 
two operators approach each other in the integral, i.e. as a: —)■ 0. This is a general feature 
in the evaluation of long-distance contributions to physical processes. In the evaluation of 
the rare kaon decay amplitude (6) Oi is the vector current and O 2 is the effective AS* = 1 
Hamiltonian density and we start by briefly recalling the normalization of these operators 
before considering the contact terms arising as they approach each other. 

A. Renormalization of Hw and 

We assumed in Sec. IIIB that the current is the conserved one given by Noether’s 
theorem applied to the chosen QCD action being used. It therefore satishes the following 
vector Ward identity: 

d^(V^(x)0(xi,... ,Xn)) = 0 (33) 

where 0(xi, ..., Xn) is a multi-local operator with all the points Xi,... ,Xn distinct from x. 
On the lattice, this identity is exactly satished and the derivative becomes a backward finite 
difference operator. This exact conservation means that the vector current does not require 
any renormalization as the continuum limit is taken. 

The Wilson coefficients Ci ^2 in the weak Hamiltonian Hw defined in (7) are currently 
known at NLO in the MS scheme [10]. Since renormalization conditions based directly on 
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Ju 



Figure 6. A potentially quadratically divergent insertion into the S and E classes of diagram. Jy 
represents the conserved electromagnetic current and Jj^ij is the local vector current Ujj^Ui or 
Cj'jf.Ci from Qi^ 2 - 


dimensional regularization, such as the MS scheme, are purely perturbative we envisage 
following the standard practice of renormalizing the bare lattice operators and 

Q2 ~ Q2 non-perturbatively into a scheme such as RI-MOM or RI-SMOM [25-27] and then 
to use continuum perturbation theory to match these renormalized operators into the MS 
scheme. 

The use of a lattice formulation with good chiral symmetry, such as domain wall fermions, 
prevents mixing with dimension-6 operators which transform under different representations 
of the chiral flavor group SU{4)i x S'f/(4)R. Within this formulation we can follow the 
renormalization procedure described in detail in [21, 22] in the evaluation of the K^-Ks 
mass difference. In that case the effective weak Hamiltonian is a simple extension of (7), 

= % E V,iVi.{CiQT' + C 2 QV') (34) 

^ q,q'=u,c 

where the operators are generalizations of those in Eq. (8) 

QV' = and Qf' = . (35) 

Since the components with q ^ q' do not contribute to the matrix elements for K —)• 
decays, one is able to rewrite Hw in Eq. (34) in the form given in Eq. (7). 

B. Additional divergences as H\y{x) approaches J(0) 

In diagrams of the “loop” class in topologies S and E [cf. Figs. 9 and 10), there are 
insertions of the form illustrated in Fig. 6. This has been studied in some detail in [1] and 
we briefly summarise the conclusions. The vector current J^, to which the photon couples is 
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the conserved one whereas the vector cnrrent from the weak Hamiltonian is a local one; 
the label L represents Local. By power connting the loop integral appears to be qnadratically 
divergent. This is reminiscent of the evalnation of the one-loop contribntion to the vacnnm 
polarization in QED and QCD and jnst as in those cases, electromagnetic gange invariance 
implies that there is a transversality factor of — q^g^'^ and the order of divergence 
is rednced by two to a logarithmic one. (In momentnm space with a lattice action the 
Ward identity q'"Jy = 0 becomes q'^Jy = 0, with = (2/a) sin(ag^/2).) This strnctnre 
was verihed and the divergence explicitly calcnlated in [1] in one-loop lattice pertnrbation 
theory for Wilson, clover and twisted-mass fermions. The logarithmic divergence is mass 
independent, and so cancels exactly in the GIM snbtraction between the diagrams with u 
and c-qnark loops. 

The above argnment can be extended straightforwardly to higher-order diagrams in which 
the gluons are contained within the n or c quark loop in Fig. 6. The emission of one or more 
gluons from the m or c propagators in the loop to be absorbed by a quark or gluon propagator 
which is external to the loop reduces the order of divergence, again leading to a convergent 
loop integration as Ju{x) approaches Hw- The remaining divergences are those which are 
associated with the renormalization of Hw- 

We have seen that as a result of gauge invariance and the GIM mechanism in the four- 
flavor theory there are no additional UV divergences in J d^x (7r| T[J(0) Hy/i^x)] \K) coming 
from the short distance region x ~ 0. In the three-flavor theory, gauge invariance still 
protects the correlation function from quadratic divergences, but then there remains a loga¬ 
rithmic term which can be removed using non-perturbative renormalization techniques [25]. 


VI. CONCLUSIONS 

Precision flavor physics will continue to be a central tool in searches for new physics and 
in guiding and constraining the construction of theories beyond the Standard Model. Lat¬ 
tice QCD simulations play an important role in quantifying the non-perturbative hadronic 
effects in weak processes. We must therefore continue to both improve the precision of the 
determination of standard quantities (such as leptonic decay constants, semileptonic form 
factors, neutral meson mixing amplitudes etc.) and to extend the range of physical quantities 
which become amenable to lattice simulations. In this paper we propose a procedure for the 
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evaluation of the long-distance effects in the rare kaon decay amplitudes K —)■ . These 

effects represent a significant (and unknown) fraction of the amplitudes. In a companion 
paper [2] we discuss the prospects for the evaluation of long distance contributions to the 
rare decays K —>■ ttvV which will soon be measured by the NA-62 experiment at CERN and 
the KOTO experiment at J-PARC. These decays are dominated by short-distance contri¬ 
butions, but given that they will soon be measured, it is interesting also to determine the 
long-distance effects which are expected to be of the order of a few percent for decays. 

In the previous sections we have explained how the technical issues needed to perform 
the lattice simulations can be resolved. Unphysical terms which grow exponentially with 
the range of the time integration, generally present when evaluating long-distance effects 
containing intermediate states with energies which are less than those of the external states, 
were shown in Sec. Ill C to be absent or small. They could potentially arise from the presence 
of intermediate states consisting of one, two or three pions and we considered each of these 
cases in turn. Similarly, the corresponding finite-volume corrections are small provided that 
the invariant mass of the lepton-pair is smaller than 2Mn-. Ultraviolet effects were discussed 
in Sec. V. We envisage using the lattice conserved electromagnetic vector current so no 
renormalization of this operator is required. In addition to the now standard renormaliza¬ 
tion of the weak Hamiltonian Hy/, we need to consider the possible additional ultraviolet 
divergences which may arise when and Hy/ approach each other. Electromagnetic gauge 
invariance implies that no quadratic divergence is present [1] and in the four-flavor theory 
the remaining logarithmic divergence is cancelled by the GIM mechanism. 

We conclude that it is now feasible to begin studying rare kaon decays K —)■ 
in lattice simulations and in particular to computing the long-distance contributions. The 
next step is to figure out how to practically implement the numerical calculation of these 
amplitudes. The main challenge resides in the calculation of the 4-point function in Eq. (18). 
One important problem is the evaluation of diagrams containing closed loops which require 
the knowledge of quark propagators from a point to itself for every lattice site. Actual 
numerical calculations have been underway for the past two years [28] and we are preparing 
papers describing exploratory results. Within the next five years we would hope that the 
hadronic effects in these decays would be quantified with good precision, thus motivating 
the extension of the experimental studies of A' —?• ttvV decays to include also K —)■ 
decays. 
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Appendix A: Feynman diagrams contributing to the rare kaon decay correlator 




IFdis. 


Figure 7. “Wing” class of diagram contributing to the rare kaon decay correlator (18). The dia¬ 
grammatic conventions are the same as those in Fig. 3. 
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Figure 8. “Connected” class of diagram contributing to the rare kaon decay correlator (18). The 
diagrammatic conventions are the same as those in Fig. 3. The £ quark is an up or down quark 
depending on the charge of the initial and final states. 
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u, c 



Figure 9. “Saucer” class of diagram contributing to the rare kaon decay correlator (18). The 
diagrammatic conventions are the same as those in Fig. 3. 


23 















Figure 10. “Eye” class of diagram contributing to the rare kaon decay correlator (18). The dia¬ 
grammatic conventions are the same as those in Fig. 3. 



Figure 11. “tt^ ” class of diagram contributing to the rare kaon decay correlator (18). These diagrams 
only contribute to the decay of the neutral meson K^. The diagrammatic conventions are the same 
as those in Fig. 3. 
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